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Abstract 

We perform a detailed study of the critical behavior of the mean field diluted Ising 
ferromagnet by analytical and numerical tools. We obtain self-averaging for the magne- 
tization and write down an expansion for the free energy close to the critical line. The 
scaling of the magnetization is also rigorously obtained and compared with extensive 
Monte Carlo simulations. We explain the transition from an ergodic region to a non triv- 
ial phase by commutativity breaking of the infinite volume limit and a suitable vanishing 
field. We find full agreement among theory, simulations and previous results. 

1 Introduction 

In the past few years a match among the study of systems defined on lattices by means of 
statistical mechanics [37] and the study of networks by means of graph theory [T3] gave origin 
to very interesting models as the small world magnets [35j or the scale free networks [13] . 
However a complete analysis starting from the simple fully connected mean field Ising model 
[6] [l9] up to these recent complex models [31] is still not complete (even thought several 
important steps have been obtained, examples being [29] [E] [21] ) and important models have 
not yet been taken into account. Among these a certain role is played by the mean field diluted 
Ising model: an Erdos-Renyi networks [5D] which has spins as nodes and their interactions as 
links. The interactions are encoded in a matrix connecting pairs of spins which, when the 
connectivity allows the link to be present, shares the same value for all the couples. 
Despite their easy formalization [35j , diluted ferromagnets are poorly investigated by rigorous 
tools [22] . Inspired by a recent work on these systems in which the authors presented a detailed 
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analysis of the ergodic region and the zero temperature hne [18] we extend recent techniques 
developed in a series of papers [3] [5] [2] [10] to this model with the aim of analyzing its critical 
behavior. We systematically develop the interpolating cavity field method 25J and use it to 
sketch the derivation of a free energy expansion: the higher the order of the expansion, the 
deeper we could go beyond the ergodic region. Within this framework we perform a detailed 
analysis of the scaling of magnetization (and susceptibility) at the critical line. The critical 
exponents turn out to be the classical ones. At the end we perform extensive Monte Carlo 
(MC) simulations for different graph sizes and bond concentrations and we compare results 
with theory. Indeed, also numerically, we provide evidence that the universality class of the 
diluted Ising model is independent of dilution. In fact the critical exponents we measured 
are consistent with those pertaining to the Curie- Weiss model, in agreement with analytical 
results. The critical line is also well reproduced. 

The paper is organized as follows: In Section ([2|) we describe the model, in Section ([3]) we 
introduce the cavity field technique, which constitutes the framework we are going to use in 
Section ^ to investigate the free energy of the system at general values of temperature and 
dilution. Section ([5]) deals with the criticality of the model; there we find the critical line 
and the critical behavior of the main order parameter, i.e. magnetization, we provide its self- 
averaging and work out a picture by which we explain the breaking of the ergodicity. Section 
([6]) is devoted to numerical investigations, especially focused on criticality. Finally, Section ^ 
is left for outlook and conclusions. 

2 Model and notations 

Given N points and families of i-i-d random variables uniformly distributed on these 

points, the (random) Hamiltonian of the diluted Curie- Weiss model is defined on Ising A^-spin 
configurations a = (<ti, . . . , cttv) through 



where is a Poisson random variable with mean C and a > 1/2 is the connectivity. The 
expectation with respect to all the (quenched) random variables defined so far will be denoted 
by E, while the Gibbs expectation at inverse temperature f3 with respect to this Hamiltonian 
will be denoted by il, and depends clearly on a and (3. We also define (•) — Er2(-). The 
pressure, i.e. minus /3 times the free energy, is by definition 




(1) 



A^ia) = ^Eln ZN{f3) = ^-Eln^ cxp(-/3Fjv(a, a)) 
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where we implicitly introduced the partition function Zn((3) too. When we omit the depen- 
dence on N we mean to have taken the thermodynamic limit which we assume to exist for all 
the observables we deal with, in particular for the free energy [52] [H] [IH] (however we will look 
for a more firm ground on this point by numerical investigation in sec.®)- The quantities 
encoding the thermodynamic properties of the model are the overlaps, which are defined on 
several configurations {replicas) cr^^\ . . . , fr*^"^ by 

1 (1) («) 

i=l 

Particular attention must be payed on gi = m = cTi which is called magnetization. 

When dealing with several replicas, the Gibbs measure is simply the product measure, with the 
same realization of the quenched variables, but the expectation E destroys the factorization. 
Sometimes for the sake of simplicity we will call 9 — tanh(/3). 



3 Interpolating with the cavity field 

In this section at first we introduce the cavity field technique on the line of [5j by expressing 
the Hamiltonian of a system made of iV + 1 spins through the Hamiltonian of N spins by 
scaling the connectivity degree a and neglecting vanishing terms in N as follows 



u=l u=l 

such that we can use the more compact expression 

HN+i{a) ^ HN{a) + HN{a)aN+i 



P2a 



(2) 



(3) 



with 



N N^oo 

-a — > a. 



P2a 



N+1 



(4) 



So we see that we can express the Hamiltonian of iV + 1 particles via the one of N particles, 

paying two prices: the first is a rescaling in the connectivity (vanishing in the thermodynamic 
limit), the second being an added term, which will be encoded, at the level of the thermody- 
namics, by a suitably cavity function as follows: let us introduce an interpolating parameter 
t G [0, 1] and the cavity function ^'jv('5,t) given by 



*(a,/3;t) 



lim ^'jv(Q!,/3;t) 



(5) 



hm E 



In 



lim ] 

N^oo 



In 



^JV,t(a,/3) 
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The three terms appearing in the decomposition ^ give rise to the structure of the following 
theorem: 

Theorem 1 In the N ^ oo limit, the free energy per spin is allowed to assume the following 
representation 

A(a,/3) =ln2-a^^^^ + ^(a,/3;t= 1) (6) 
oa 

Proof 

Consider the A'^ + 1 spin partition function Zpf^i(a, f3) and let us split it as suggested by eq. 

m 

ZN+i{a,(3) = Yl e-P"-^^^"^^ e-''^"(")"''^"(")^«+i (7) 

{c^Af + l} {c^Af + l} 

{fN + l} 

where the factor two appears because of the sum over the hidden ctat+i variable. Defining a 
perturbed Boltzmann state a) (and its replica product = tj x ... x u)) as 



E{o-„}5(^)e 



Z^{(Tjv} i 

where the tilde takes into account the shift in the connectivity a — *■ a and multiplying and 
dividing the r.h.s. of eq.® by ZN{a,l3), we obtain 

ZN+i{a,P) = 2ZN{a,l3)io{e^^^='i). (8) 

Taking now the logarithm of both sides of eq. ([8]) , applying the average E and subtracting the 
quantity [In Zpf+iio-, P)], we get 

E[lnZA,+i(a,/3)] - E[ln Za,+i(5, ^)] ^ In2 + E[ln /^^"l^^ l + N^a, l3;t ^ I) (9) 



ZN+iia,f3) 



in the large N limit the l.h.s. of eq.® becomes 



E[lnZjv+i(a,/?)] -E[lnZAr+i(a,/3)] = (10) 
[a- a)—W.[hiZN+i{a,(i)\ = a ^ ^ ^ —[hiZN+i{a, = a — 

then by considering the thermodynamic limit the thesis follows. (Actually we still do not 

have a complete proof of the existence of the thermodynamic limit but we will provide strong 
numerical evidences in Section [6]) □ 



Hence, we can express the free energy via the energy and the cavity function. While it is 
well known how to deal with the energy [18j . the same can not be stated for the cavity function, 
and we want to develop its expansion via suitably chosen overlap monomials in a spirit close 
to the stochastic stability [3][IS][33], such that, at the end, we will not have the analytical 
solution for the free energy in the whole (a, /3) plane, but we will manage its expansion close 
(immediately below) to the critical line. To see how the machinery works, let us start by giving 
some definitions and proving some simple theorems: 

Definition 1 We define the t-dependent Boltzmann state ujt as 



where Zjv, (a/?) extends the classical partition function in the same spirit of the numerator of 



As we will often deal with several overlap monomials let us divide them among two big cate- 



Definition 2 We can split the class of monomials of the order parameters in two families: 

• We define filled or equivalently stochastically stable all the overlap monomials built by 
an even number of the same replicas (i.e. qi2, m? , (7i2'J'34'?i234/'- 

• We define fillable or equivalently saturable all the overlap monomials which are not 
stochastically stable (i.e. qi2, m, gi2'?34j 

We are going to show three theorems that will play a guiding role for our expansion: as this 
approach has been deeply developed in similar contexts (as fully connected Ising model [B] 
or fully connected spin glasses [5] or diluted spin glasses [8] , which are the boundaries of the 
model of this paper) we will not show all the details of the proof, but wc sketch them as they 
are really intuitive. The interested reader can deepen this point by looking at the original 
works. 

Theorem 2 For large N , setting t — 1 we have 



such that in the thermodynamic limit, if t = 1, the Boltzmann average of a fillable multi- 
overlap monomial turns out to be the Boltzmann average of the corresponding filled multi- 
overlap monomial. 




(11) 



e<?.(Qip. 



gories: 



'^Ar,f((Tii(Ti2...(Ti„) u)N+i{<Ji^ai^...ai^<j'^^i) + 0(— ) 



(12) 
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Theorem 3 Let Q2n be a fillahle monomial of the overlaps (this means that there exists a 
multi- overlap q2n such that q2nQ2n *s filled). We have 

lim lim((52n)t = {q2nQ2n) (13) 

(examples: for N ^ oo we get (mi)* (ml), (912)* {^12)1 (912934)* (91293491234) J 

Theorem 4 In the N ^ 00 limit the averages (•) of the filled monomials are t-independent 
in P average. 

Proof 

In this sketch we are going to show how to get Theorem ^ in some details; It automaticaUy 
has as a corollary Theorem (jSj which ultimately gives, as a simple consequence when applied 
on filled monomials, Theorem Q. 

Let us assume for a generic overlap correlation function Q, of s replicas, the following repre- 
sentation 

Q-flEfl<r^(W» 

a— 1 l—l 

where a labels the replicas, the internal product takes into account the spins (labeled by 
/) which contribute to the a-part of the overlap qa.a' and runs to the number of time that 
the replica a appears in Q, the external product takes into account all the contributions of 
the internal one and the / factor fixes the constraints among different replicas in Q; so, for 
example, Q — (712923 can be decomposed in this form noting that s = 3, = n"^ = = 2, 
/ = A'^~^(5ji j3(5j2 j3, where the 6 functions fixes the links between replicas 1,2 —> qi^2 and 
2, 3 — !■ (72. 3- The averaged overlap correlation function is 

(g),=E5]/(K})n'^*(n<?)- 

i'^ a—1 l — l 

Now if Q is a fillable polynomial, and we evaluate it at t = 1, let us decompose it, using the 
factorization of the uj state on different replica, as 

ij>^ii> a=l 1=1 b=u 1=1 

where u stands for the number of the unfilled replicas inside the expression of Q. So we split 
the measure VL into two different subset oja and uji,: in this way the replica belonging to the h 
subset are always in even number, while the ones in the a subset are always odds. Applying 
the gauge erf — > crf(T^^]^,Vi G (l,iV) the even measure is unaffected by this transformation 



6 



(f^^+i — 1) while the odd measure takes a (Jn-^i inside the Boltzmann measure. 

u n°' s 

(Q) = E ^({*"}' n ^(^^^+1 n <r) n n 4) 

a=l /=1 b=u Z=l 

At the end we can replace in the last expression the subindex TV + 1 of (Tat+i by k for any 
k ^ {z"} and multiply by one as 1 = 7V-1 Up to orders 0{1/N), which go to zero m 

the thermodynamic limit, we have the proof. □ 

It is now immediate to understand that the effect of Theorem ([5]) on a fillable overlap monomial 
is to multiply it by its missing part to be filled (Theorem[3]), while it has no effect if the overlap 
monomial is already filled (Theorem |3|) because of the Ising spins (i.e. (jj^^i = 1 Vrt G N). 

Now the plan is as follows: We calculate the f-streaming of the function in order to 
derive it and then integrate it back once we have been able to express it as an expansion in 
power series of t with stochastically stable overlaps as coefficients. At the end we free the 
perturbed Boltzmann measure by setting t — 1 and in the thermodynamic limit we will have 
the expansion holding with the correct statistical mechanics weight. 

^5%M ^ |E[ln^(e^Sri!'-S] (14) 
dt dt ^ ^ ^' ^ ' 

= 25E[ln(:;;(e'^^-''i''^'-+'^'''o)] -2aE[lnc:;(e'^^-i*'^"0] = 25E[ln J)t(e'''"'o )] 
using now the equality e^°'^o = cosh/3 + (Tig sinh/3, we can write the r.h.s. of eg. p^ as 

d t) 

-^f^ = 25E[lnwt(cosh/? + cr,o sinh/3)] = 2(51og cosh/3 - 2(5E[ln(l + ^^((CTije)]. 

We can expand the function log(l + uitS) in powers of 6, obtaining 

= 25 In cosh/3 - 25 V (15) 
ot ^-^ n 

n— 1 

We learn by looking at eq. (fT5|) that the derivative of the cavity function is built by non- 
stochastically stable overlap monomials, and their averages do depend on t making their t- 
integration non trivial (we stress that all the fillable terms are zero when evaluated at t = 
due to the gauge invariance of the model). We can escape this constraint by iterating them 
again and again (and then integrating them back too) because their derivative, systematically, 
will develop stochastically stable terms, which turn out to be independent by the interpolating 
parameter and their integration is straightforwardly polynomial. To this task we introduce 
the following 
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Proposition 1 Let be a function of s replicas. Then the following streaming equation holds 

a=l 



dt 



(16) 



+ 2a^2E(F.«),, - sJ2{F,alat^\, + '±±11 (F^at:' at:%,] 



a<b 

l,s 



a=l 



l,s 



a<6<c 



o<6 



+ 



s{s + 1)m-^, „ I S{s + 1){S + 2) ^s+3\ 1 



2! 



a=l 



where we neglected terms 0{6^). 



Proof 

The proof works by direct calculation: 



d{Fs)t,a 

dt 



2<5E 



2dE 



25E 



= 2aE 



(17) 



fit(e 



f}t(F,n^=i(cosh/3 + <7? sinh/3)) 



^^t(n^=i(cosh/3 + c7« sinh/J)) 

f]t(F,n«^i(i + <0)) 



2a{F.s) 



t.a 



(l+d;t«)0)^ 



- 2a(F«)t,5 = 



Now noting that 



1 



a=l 



a<6 



a<6<c 



we obtain 

at 



(i + wt^)^ 



25E 



2! 



3! 



a=l o<6 o<6<c 



(18) 



(19) 



from which our thesis follows □. 
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4 Free energy analysis 



Now that we exploited the machinery we can start applying it to the free energy. Let us at 
first work out its streaming with respect to the plan (a,/?): 



dA[a,P) (H) 1/1 



dp N iV , , , 

{<t} f=l 



= ^^kTr{k - l,aN)E[Lj{a,^ajJk] = a ^ 7r(fc - 1, aiV)E 

k=l k=l 

L a;(cosh/? + (Tij.crj^ sinh/3) J 11 + uj{ai^aj^)9 



by which we get (and with similar calculations for daA{aT0) that we omit for the sake of 
simplicity) : 

^ ae-af2{-ini-0')e-~'{ql.,J (21) 

n— 1 

. ,„».t„-f Hrr(,t,..) (22) 

71—1 

Now remembering Theorem ([l} and assuming critical behavior (that we will verify a fortiori 

in sec. ([5])) we move for a different formulation of the free energy by considering the cavity 
function as the integral of its derivative. In a nutshell the idea is as follows: Due to the second 
order nature of the phase transition for this model (i.e. criticality that so far is assumed) we 
can expand the free energy in terms of the whole series of order parameters. Of course it is 
impossible to manage all these infinite overlap correlation functions to get a full solution of 
the model in the whole (a, (3) plane but it is possible to show by rigorous bounds that close 
to the critical line (that we are going to find soon) higher order overlaps scale with higher 
order critical exponents so we are allowed to neglect higher orders close to this line and we 
can investigate deeply criticality, which is the topic of the paper. 
To this task let us expand the cavity functions as 

f(a,/3,i) = —dt' (23) 

= 25ilogcosh/3 + /3 / {m)t',&dt' -]-pe I {qi2)t',&dt' + 0{e^) 
Jo ^ Jo 

where /3 = 2olQ ^ (3' = 2a9 for — > oo. Now using the streaming equation as dictated by 
Proposition ([1]) we can write the overlaps appearing in the expression of 5* as polynomials of 
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higher order filled overlaps so to obtain a straightforward polynomial back-integration for the 
5* as they no longer will depend on the interpolating parameter thanks to Theorem ([4]). 
For the sake of simplicity the a-dependence of the overlaps will be omitted keeping in mind 
that our results are all taken in the thermodynamic limit and so we can quietly exchange a 
with a in these passages. The first equation we deal with is: 

= (3[{m') - (mima)*] (24) 

dt 

where (miTO2) is not filled and so we have to go further in the procedure and derive it in order 
to obtain filled monomials: 

= 2/3[(mim2)t-(miTO2m3)t]-|-/36l[(mim2(7i2)-4(TOim2(?i3)t+3(miTO2q34)t]- (25) 

In this expression we stress the presence of the filled overlap (mim2qi2) and of (771^7712)4 which 
can be saturated in just one derivation. Wishing to have an expansion for {m)t up to the third 
order in 9, it is easy to check that the saturation of the other overlaps in the last derivative 
would carry terms of higher order and so we can stop the procedure at the next step 

^^""Ir^^* ^[(r77?m2)] + /3[unfilled terms] + 0(6*^) (26) 
from which integrating back in t 

{mlm2)t = P[{mlml)]t (27) 
inserting now this result in the expression (j25p and integrating again in t we find 

(77ii7772)t = pe{mim2qi2)t + P'^{m\ml)t^ (28) 
and coming back to {m)t we get 

{m)t = (3{m^)t - ^(mi7772gi2)i' - ^{mlml)t'' (29) 
which is the attempted result. Let us move our attention to ((712)4, analogously we can write 
= 2/3[(7r7igi2)4 - (7713912)4] + /3e[(g?2> - 4(qi2gi3)t + 3(912934)4] (30) 

dt 

and consequently obtain 

(912)4 = P0{ql2)t + P\m^m2qi2)t' + 0(0^). (31) 
With the two expansion above, in the iV — > 00 limit, putting t = 1 we have 
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*(a,/3,i = 1) = 2alncosh/3+ ^(m^) - ^{mlml) - ^{ql^) - ^(mimagia) + 0{e^) 

(32) 

At this point we have aU the ingredients to write down the polynomial expansion for the free 
energy function as stated in the next: 



Proposition 2 A general expansion via stochastically stable terms for the free energy of the 
diluted Ising model can he written as 

3' 

A{a,(3) = ln2 + alncosh/3+y (/3'-l)(m?)+ (33) 



It is immediate to check that the above expression, in the ergodic region where the averages 
of all the order parameters vanish, reduces to the well known high-temperature (or high con- 
nectivity) solution [18] (i.e. A{a,3) = In 2 + a log cosh /3) . 

Of course we are neglecting 6^ higher order terms because we are interested in an expansion 
holding close to the critical line, but we are not allowed to truncate the series for a general 
point in the phase space far beyond the ergodic region. 

5 Critical behavior 

Now we want to analyze the critical behavior of the model: we find the critical line where the 
ergodicity breaks, we obtain the critical exponent of the magnetization and the susceptibility 
and at the end we show that within our framework the lacking of ergodicity can be explained 
as the breaking of commutativity of the infinite volume limit against our cavity field, thought 
of as a properly chosen field, vanishing in the thermodynamic limit too, accordingly with the 
standard prescription of statistical mechanics [T] . 

5.1 Critical line 

Let us firstly define the rescaled magnetization as ^^r = \/NmN. By applying the gauge 
transformation at ~^ aiCFN+i in the expression for the quenched average of the magnetization 
(eq. (HH])) and multiplying it times N so to switch to setting t = 1 and sending N ^ oo 
we obtain 

(e?> = ^^^^(Ci^mim^) + ^^|^(a6gi2) + O(^) (34) 
by which we see (again remembering criticality and so forgetting higher order terms) that 
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10| , , , , , 

g . I percolation line 

' critical line 



Figure 1: Phase diagram: below ac = 0.5 there is no giant component in the Erdos-Renyi 
graph, ac defines the percolation threshold. Above at left of the critical line the system behaves 
ergodically, conversely on the right ergodicity is broken and the system displays magnetization. 



the only possible divergence of the (centered and rescaled) fluctuations of the magnetization 
happens at the value P' — 1 which gives 2a6 = 1 as the critical line, in perfect agreement with 
[18]. The same critical line can be found easier simply by looking at the expression (j33p as 

follows: remembering that in the ergodic phase the minimum of the free energy corresponds 
to a zero order parameter (i.e.-\/ {m'^) = 0), this implies the coefficient of the second order 
a{P') = ^{P' — 1) to be positive. Anyway immediately below the critical line values of the 
magnetization different from zero must be allowed (by definition otherwise we were not crossing 
a critical line) and this can be possible if and only if a(/3') < 0. Consequently (and using once 
more the second order nature of the transition) on the critical line we must have a{(3') = 
and this gives again 2a6 = 1. 

5.2 Critical exponents and bounds 

Now let us move to the critical exponents: 

Critical exponents are needed to characterize singularities of the theory at the critical line and, 
for us, these indexes are the ones related to the magnetization (rn) and to the susceptibility 

(x). 

We define r — (2Q;tanh/3 — 1) and we write {m{T)) ^ Go ■ and {x{t)) ~ Go • t'*', where the 
symbol ^ has the meaning that the term at the second member is the dominant but there are 
corrections of higher order. 
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Remembering the expansion of the squared magnetization that we rewrite for completeness 
(™^) = Tri^ — Tt("^i™2) + 777^~~TT ('^1^2912) + 0{ J ) (35) 



and considering that using the same gauge transformation ai —t aia^^i on (cq. (j3ip l we have 
for the two rephca overlap the following representation 

{QI2) = - (^,f^^) (^i^2gi2) + 0{e') (36) 

we can arrive by simple algebraic calculations to write down the free energy, close to the critical 
line of course, depending only by the two parameters {m^) and 

Aia, P) = ln2 + alncosh/3 + ^ (/?' - 1) (m?) - ^ - 1^ (q^^) + 0(^6) (37) 



By a comparison of the formula obtained by deriving A(a, (3) as expressed by eq. p7p and the 
expression we have previously found (eq. ([55)) ) that we report for the sake of readability. 



dA{a,l3) 
da 



n=l 



= In cosh /? - V —0"(ql 

^-^ n 



(38) 



it is immediate to see that we have 



(39) 



If we put ourselves close to the value P' — I and make a change of variable r = /?' — 1 with 
da = 26dr we get 

^ r/5',«/ ^^/ 2x1 ^ ^ 



da 



by which we easily obtain 



2dT 



Ot d{ml) 



2 dr 



e{mi) 



d{ra\) dr 
(m?) T 



(40) 



(41) 



Therefore we get that the critical exponent for the magnetization, 5 — 1/2, which turns out 
to be the same as in the fully connected counterpart [B] [19] , in agreement with the disordered 
extension of this model [10]. 

Again, by simple direct calculations, once we get the critical exponent for the magnetization 
it is straightforward to show that the susceptibility (x) [1] obeys 



(42) 
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close to the critical line, by which we find its critical exponent to be once again in agreement 
with the classical fully connected counterpart yij . 

Now we want to show some wrong results which a naive calculation would suggest so to 
emphasize the importance of the bounds relating different monomials that we are going to 
discuss immediately later [9] |34j . Then in the next subsection, we explain what is the physics 
behind this picture by providing a mechanism for the breaking of the ergodicity. 
The point on which we focus is the following: if we wish to perform the same procedure we 
performed on (m^), applying blindly saturability below to the first critical line, to the 2 replica 
overlap ((712) would gain 

^) -T^^{f3'e-l) (43) 

identifying 9^^ — \/{2a)^/'^ as another critical temperature, or better, the critical temperature 
typical of ((712) ■ In the same way we could find 0c„ for every „ obtaining 

= l/(2a)i/" 

such that, at the end, we obtain a scenario with several transition lines, one for every order 
parameter. 

This is not a possible scenario, as generally explained for instance in [7] and as dictated, in 
this model, by the following 

Proposition 3 As soon as the first order parameter (the magnetization) starts taking values 
different from zero, the same happens to all the other order parameters 

{ql J = EAK{<J^,a,,)] > {EkEM'^^,<^^2)])" = (™'>" Vn (44) 
We omit the proof details as they are a simple application of Jensen inequality. 

5.3 Saturability breaking 

So far we showed that it is not possible to have several transition lines, one for every order 
parameter. Now we want to understand why there is just one critical line by applying the 
theory developed in [9] to this model. 
Starting from Theorem ([T|) that we recall for simplicity 

A(a,/?) =ln2-a^^^^^ + *(a,/3,i= 1) (45) 
oa 

we want to show the phase transition expressed by the non-commutativity among the thermo- 
dynamic limit and the vanishing perturbation. 
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Again for simplicity we report the expansion of A{a,(3) we have previously built 

Aia, (3) = ln2 + alncosh/3 + ^ (/?' - 1) (m?) -^(^-l] (qf^) + 0(^6) (46) 



4 ' " " 48a V 2a 
that we obtained by considering the cavity function as the integral of its t-derivative 



= J* ^^^j^dt' = 2a\ncosh /3dt' - 25^^ ^T"^" f 



[qi,...,n)tdt' (47) 



and performing, via the streaming equation, a saturating procedure of consecutive i-derivatives 
upon the overlaps in order to express them as functions of higher order filled terms. Then the 
only thing we had to do was sending N to infinity, carrying out of the integral the overlaps and 
then putting < = 1 to evaluate 5'(Q;,/3,t = 1). The result of this procedure brings to equation 
(gSl). 

But what if we exchanged the limit order by sending t ^ 1 first and taking the thermodynamic 
limit after? 

It is easy to note that all the overlaps appearing in eq. (|47|l are fiUable such that we can avoid 
the saturation procedure simply by setting t = 1 first and then sending N to infinity. In this 
way, thanks to theorem ([3]), each fiUable overlap is transformed in a filled ^-independent one 
and this kills all the correlations among different replicas and allows us to write 

A{a,f3) =ln2 + alncosh (3- Y ^^^e^'iql „)=ln2 + a^^^^ (48) 

n—l 

being clearly 

lim lim4'(a,/3,t) ^ lim lim [ ^Ijfi^iMldt' (49) 

{qi....,n)'tdt'] 



lim limf / 2(31ogcosh/3dt' - 2a V '—6" / 



2alncosh/3-2a^^^— „) ^ 2a 



dc 



In particular, retaining just the first two terms of the expansions, we show the difference 
between the two results. 

• limf^i limjv-,00 

B' 

A{a,(3) = ln2 + alncosh/3 + J (/3' - 1) (to?) + O{0^) (50) 

• limAT^oo hnif^i 

B' 

A(a,/3) = ln2 + alncosh/3 + ^(m?) +0(6*3) (-gj^^ 
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We immediately recognize that eq. ([501 is the correct expression holding also below the critical 
line. When the system lives in the ergodic region all the order parameters are zero and it 
reduces to Q;lncosh(/3) which is the ergodic solution, if we cross the critical line the formula 
takes into account the phase transition encoded in the coefhcient of the second order and gives 
the correct expression immediately below. 

It is also straightforward to recognize as the ergodic solution eq. ([5T|) which can be the correct 
one only above the critical line [TH] . 

At the end we saw that there exists only and only one critical line for all the order parameters. 
We saw that this line can be depicted as the breaking of commutativity among the infinite 
limit operation and the setting of i = 1 (relaxing the Boltzmann factor from the interpolant 
avoiding the trivial way t = 0). Hence, which can be the genesis of the transition for all 
the other order parameters? The correlations among them and the magnetization, as clearly 
put in evidence in the free energy expansion (see ea. (j33p V And lastly, which is the origin of 
these correlation? The saturability property that the model exhibits as stated by Theorem ([3]). 
In fact, at the critical line (which is the last line, from above, where gauge invariance is a 
symmetry of the Boltzmann state too, thanks to the continuity of the transition) saturability 
easily shows that 

lim {miqi2)t = ("11^2912), (52) 

explaining the birth of the correlations among the various order parameters (we reported just 
the first two as an example). 

5.4 Self-averaging properties 

We have previously shown how filled overlaps become asymptotically i-independent when N 
grows to infinity. Starting from this we can find identities stating the self-averaging of the 
order parameters as (m?). 

In particular we are going to take these overlaps and calculate their derivative with respect to 
t; then by applying the gauge transformation and setting i = 1,7V ^ 00 we can write down 
the self- averaging relations. 

= dt{ml)t = P'[{ml)t - {mlm^),] [{mi) - {mlml)] = m{m^) - ^\m^)]. (53) 

Coherently with [18j we find standard self-averaging for the magnetization. 
More interesting is the situation concerning the overlap but we actually lack a complete math- 
ematical control. 

By applying the same trick as before we can equate to zero (in the large N limit) the t- 
derivative of the squared overlap (which is stochastically stable), consequently we get two 
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terms. 



- 9t(g?2>t = 2(3'[{m,qi^)t ~ (m3g?2>t] + /^'^^K'zD* - 4(g?29i3)t + 3(q?2'?34)t] (54) 
^ [(™?9?2> - (™3'Z?2>] = [(912) - 4(9?2'Z?3> + 3(g?2'?3'4>] = (55) 

A-priori we can not assume factorization of the series so to put to zero each term separately (as 
we did in eq. (|55p ). however, close to the critical line surely the second term is an higher order 
and we can reasonably set to zero the first. Furthermore as the second term at the r.h.s. has 
a pre-factor cx a^^ differing from the first term it would be difficult to imagine the opposite. 
It is in fact very natural to assume that the two terms can be set to zero separately in the whole 
(a, P) plane and this is very interesting because the second term is a very well known relation in 
the field of spin glasses [IS] [H] [3] [13] |5] suggesting a common structure among different kinds 
of disordered systems, the only sharing feature among diluted fcrromagnets and spin-glasses 
being some kind of disorder (topological in the former, frustrating in the latter). 
We are not going to deepen this point as it is under investigation in [IT where the same set 
of relations (and more) is found and discussed. 

6 Numerics 

In this section we analyze, from the numerical point of view, the ferromagnetic system previ- 
ously introduced by performing extensive Monte Carlo simulations with the Metropolis algo- 
rithm [32]. 

The Erdos-Rcnyi random graph is constructed by taking N sites and introducing a bond 
between each pair of sites with probability p ~ a/{N — 1), in such a way that the average 
connectivity per node is just a. Clearly, when p = 1 the complete graph is recovered. 
The simplest version of the diluted Curie- Weiss Hamiltonian has a Poisson variable per bond 
as — — X^ij X^r^o" '^iu'^ji^j ^^'^ it is the easiest approach when dealing with numerics. 
For the analytical investigation we choose a slightly changed version (see eq.([T])): each link 
gets a bond with probability close to a/N for large N; the probabilities of getting two, three 
bonds scale as therefore negligible in the thermodynamic limit. 

Working with directed links (as we do in the analytical framework) the probability of having 
a bond on any undirected link is twice the probability for directed link (i.e. 2a /N). Hence, 
for large N, each site has average connectivity 2a. Finally in this way we allow self-loop but 
they add just cr-independent constant to Hj^ and are irrelevant, but we take the advantage of 
dealing with an Hn which is the sum of independent identically distributed random variables, 
that is useful for analytical investigation. 
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When comparing with numerics consequently we must keep in mind that a = 2a. 
In the simulation, once the network has been diluted, we place a spin ai on each node i and 
allow it to interact with its nearest-neighbors. Once the external parameter f3 is fixed, the 
system is driven by the single-spin dynamics and it eventually relaxes to a stationary state 
characterized by well-defined properties. More precisely, after a suitable time lapse to and for 
sufficiently large systems, measurements of a (specific) physical observable x{a, a, (3) fluctuate 
around an average value only depending on the external parameters /3^^ and a. 
Moreover, for a system [a, (3) of a given finite size N the extent of such fluctuations scales 
as iV~5 with the size of the system. The estimate of the thermodynamic observables (x) is 
therefore obtained as an average over a suitable number of (uncorrelated) measurements per- 
formed when the system is reasonably close to the equilibrium regime. 

The estimate is further improved by averaging over different realizations of the same system 
{a, (3). In summary, 

1 1 

^ a;(CT(t„)) , tr, =tn + nT 

n=l 



where a{t) denotes the configuration of the magnetic system at time step t and T is the 
decorrelation parameter (i.e. the time, in units of spin flips, needed to decorrelate a given 
magnetic arrangement). 

In general, statistical errors during a MC run in a given sample result to be significantly 
smaller than those arising from the ensemble averaging (see also P?]). Figure ^ shows the 
dependence of the macroscopic observables (m) and (e) from the size of the system; values are 
obtained starting from a ferromagnetic arrangement, at the normalized inverse temperature 
/3/a — 1.67. Notice that at this temperature the system composed of = lO'* is already very 
close to the asymptotic regime. Analogous results are found for different systems (a,/3), with 
P far enough from Pc- 

In the following we focus on systems of sufficiently large size so to discard finite size effects. 
For a wide range of temperatures and dilutions, we measure the average magnetization (m) 
and energy (e), as well as the magnetic susceptibility calculated as 

Their profiles display the typical behavior expected for a ferromagnet and, consistently 
with the theory, highlight a phase transition at well defined temperatures Pdct)- 
Now, we investigate in more detail the critical behavior of the system. We collect accurate data 
of magnetization and susceptibility, for different values of a and for temperatures approaching 
the critical one. These data are used to estimate both the critical temperature and the critical 
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Figure 2: Finite size scaling for the magnetization and the internal energy (inset) for a = 10 
and ^ = 1.67. All the measurements were carried out in the stationary regime and the 
error bars represent the fluctuations about the average values. We find good indication of 
the convergence of the quantities on the size of the system and thus of the existence of the 
thermodynamic limit. 

exponents for the magnetization and susceptibility. In Fig. ([3]) we show data as a function 
of the reduced temperature r — {\p — Pd/Pc)^^ for a — 10 and a = 20. The best fit for 
observables is the power law 

(m) ^ /? > (56) 

X - r\ (57) 

We obtain estimates for /3c(q:), 5{a) and 7(a) by means of a fitting procedure. Results are 
gathered in Tab. [6] Within the errors (< 2% for (ic and < 5% for the exponents), estimates 



a 




5 


7 


10 


9.93 


0.48 


-0.97 


20 


19.92 


0.49 


-1.04 


30 


29.98 


0.48 


-1.04 


40 


39.59 


0.50 


-1.02 



Table 1: Estimates for the critical temperature and the critical exponents 5 and 7 obtained 
by a fitting procedure on data from numerical simulations concerning Ising systems of size 
N — 36000 and different dilutions (we stress that analytically we get 5 — 0.5 and 7 = — 1). 
Errors on temperatures are < 2%, while for exponents are within 5%. 

for different values of a agree and they are also consistent with the analytical results exposed 
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Figure 3: Log-log scale plot of magnetization (main figm'e) and susceptibility (inset) versus the 
reduced temperature r — (|/3 — Pd/ Pc)^^ for a ~ 10. Symbols represents data from numerical 
simulations performed on systems of size N = 36000, while hues represent the best fit. 

in Sec. ^ 

We also checked the critical line for the ergodicity breaking, again finding optimal agreement 
with the criticality investigated by means of analytical tools. 

7 Conclusions 

In this paper we developed the interpolating cavity field technique for the mean field diluted 
ferromagnet. Once the general framework has been built we used it to analyze criticality: we 
found analytically the critical line and the critical exponent of the magnetization, whose self- 
averaging is also proved. We present an argument to explain the transition from an ergodic 
phase to a broken ergodicity phase by the breaking of commutativity of two limits, volume 
and applied field, as dictated by standard statistical mechanics. We furthermore showed the 
existence of only one critical line where all the multi-overlaps start taking positive values as 
soon as the magnetization becomes different from zero. We proved this both mathematically 
by a rigorous bound and physically by a mechanism that generates strong correlations among 
magnetization and overlaps at the (unique) critical line: saturahility. 

At the end a detailed numerical analysis of the model is presented: by sharp Monte Carlo 
simulations the convergence of the energy density (and the magnetization) to its limit is 
investigated obtaining monotonicity in the system size. The critical line, as well as scaling of 
the magnetization and the susceptibility, are also investigated obtaining full agreement among 
theory and simulations. 

Future works should extend these techniques to several lateral models as the bipartite diluted 
mean field Ising models, while the need of stronger techniques to go well beyond the critical 
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line is also to be satisfied as well as their practical application to social science or biological 
networks. We plan to follow these research lines in the future. 
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